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Abstract 

We have recently solved the inverse spectral problem for one- 
parameter families of vector fields, and used this result to construct 
the formal solution of the Cauchy problem for a class of integrable non- 
linear partial differential equations in multidimensions, including the 
second heavenly equation of Plebanski and the dispersionless Kadomt- 
sev - Petviashvili (dKP) equation, arising as commutation of vector 
fields. In this paper we make use of the above theory i) to construct 
the nonlinear Riemann-Hilbert dressing for the so-called two dimen- 
sional dispersionless Toda equation (exp(ip)) tt = tpCi&i elucidating 
the spectral mechanism responsible for wave breaking; ii) we present 
the formal solution of the Cauchy problem for the wave form of it: 
(exp(ip)) tt = ip xx + ip yy ; iii) we obtain the longtime behaviour of the 
solutions of such a Cauchy problem, showing that it is essentially 
described by the longtime breaking formulae of the dKP solutions, 
confirming the expected universal character of the dKP equation as 
prototype model in the description of the gradient catastrophe of two- 
dimensional waves; iv) we finally characterize a class of spectral data 
allowing one to linearize the RH problem, corresponding to a class of 
implicit solutions of the PDE. 
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1 Introduction 



It was observed long ago pQ that the commutation of multidimensional vector 
fields can generate integrable nonlinear partial differential equations (PDEs) 
in arbitrary dimensions. Some of these equations are dispersionless (or quasi- 
classical) limits of integrable PDEs, having the dispersionless Kadomtsev 
- Petviashvili (dKP) equation u xt + u yy + {uu x ) x = [2], [3] as universal 
prototype example, they arise in various problems of Mathematical Physics 
and are intensively studied in the recent literature (see, f.i., [I] - [IE]). In 
particular, an elegant integration scheme applicable, in general, to nonlinear 
PDEs associated with Hamiltonian vector fields, was presented in [7J and a 
nonlinear d - dressing was developed in [13J. Special classes of nontrivial 
solutions were also derived (see, f.i., [12], [T5]). 

The Inverse Spectral Transform (1ST) for 1-parameter families of mul- 
tidimensional vector fields has been developed in [19] (see also [20J). This 
theory, introducing interesting novelties with respect to the classical 1ST for 
soliton equations [2U E2] , has allowed one to construct the formal solution 
of the Cauchy problems for the second heavenly equation [23] in [19] and for 
the novel system of PDEs 

U x t ~t~ ^yy ~t~ i^>U>x)x ~t~ V x U X y V y U xx 0, i^-^) 

V X t ~\~ Vyy ~\~ UV XX ~~\~ V X V X y VyV XX 0, 

in [21]. The Cauchy problem for the v = reduction of ([T]), the dKP equation, 
was also presented in [23] , while the Cauchy problem for the u = reduction 
of ([1]), an integrable system introduced in [16], was given in [25]. This 1ST 
and its associated nonlinear Riemann - Hilbert (RH) dressing turn out to be 
efficient tools to study several properties of the solution space, such as: i) 
the characterization of a distinguished class of spectral data for which the 
associated nonlinear RH problem is linearized and solved, corresponding to a 
class of implicit solutions of the PDE (as it was done for the dKP equation in 
[26] and for the Dunajski generalization [27] of the second heavenly equation 
in [28]); ii) the construction of the longtime behaviour of the solutions of the 
Cauchy problem [26] ; iii) the possibility to establish whether or not the lack 
of dispersive terms in the nonlinear PDE causes the breaking of localized 
initial profiles and, if yes, to investigate in a surprisingly explicit way the 
analytic aspects of such a wave breaking (as it was recently done for the 
(2+l)-dimensional dKP model in [26\). 
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In this paper we make use of this theory to study another distinguished 
model arising as the commutation of vector fields, the so-called 2 dimensional 
dispersionless Toda (2ddT) equation 



C1C2 



(e*) t , = 0(Ci,C2,*), (2) 



or 



^CiC 2 = (e v ) tt , V = <t>t, (3) 



also called Boyer-Finley [29] equation or SU(oo) Toda equation [30], the 
natural continuous limit of the 2 dimensional Toda lattice 



( e **i-*» _ e A>-A>-i) ; = 0n(Ci) Ca ). (4) 



The 2ddT equation was probably first derived in [31] as an exact reduction 
of the second heavenly equation; then in [32J as a distinguished example of 
an integrable system in multidimensions. Some of its integrability properties 
have been investigated in |33j and the integration method presented in [7J 
is applicable to it. Both elliptic and hyperbolic versions of (T5]) are relevant, 
describing, for instance, integrable 7i-spaces (heavens) [SDJ [SB] and integrable 
Einstein - Weyl geometries [37] -[SB], [3D]. String equations solutions [DJ of it 
are relevant in the ideal Hele-Shaw problem JS9j HDJ STJ H21 US] . 

The integrability of ([2]) follows from the fact that (|2J) is the condition of 
commutation [L\, L2] = for the following pair of one-parameter families of 
vector fields [B]: 

U = d (l + \vd t + (-Xv t + %i) \d\, 

L 2 = % + \- l vd t + (a~V - \d x , 5 

where 

v = (6) 

and A G C is the spectral parameter, implying the existence of common 
eigenfunctions of both vector fields; i.e., the existence of the Lax pair: 



(7) 
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(8) 



Equations ([7j) and (j2J) can be written in the following " Hamiltonian" form 

& 

^Ci + {^i> ^}(A,t) = 0, ^2 + {^2, V'}(A,t) = °> 

^1(2 ~^-Kl ~ {^l5^2}(A,t) = 0, 

where 

Wl = A«-% W 2 = -A-W% (9) 

and 

{f,g} { x,t)-=\{fx9t-ft9\). (10) 

If, in particular, 

& = * = — ^ C 2 = * = — g-S xjGl, (11) 
equation (j2j) becomes the following nonlinear wave equation 

0xx + ra = (e*) t , (12) 

or 

<Pxx + <P w = (e' p ) tt . (13) 

In addition, if |A| = 1 and (f) 6 R, equations ([9]) give the "real" Hamiltonian 
formulation 1 32 1 



ip x + {Hx,ij}}(fi t t) = 0, ip y + {H 2 , 4>}(e,t) = 0, 
H\ y — H 2x — {Hi, H 2 }(e,t) = 0, 

of (jTJ) and f|T2|) . for the real Hamiltonians Hi,H 2 : 

Hi = %(Hi + H 2 ) = sm9e^ - \<j) y , 
H 2 = \{H 2 -Hi) = - cosOe^ + \4> x 

and the Poisson bracket 



(14) 



(15) 



{f,9}($,t) ■= fe9t~ ftgo, (16) 

having introduced the parametrization 

\ = e - ie , 6eR. (17) 

The paper is organized as follows. In §2 we present the dressing scheme 
for equations (|2J) and (|12|) . given in terms of a vector nonlinear RH problem. 
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As for the dressing of dKP presented in since the normalization of the 
eigenfunctions turns out to depend on the unknown solution of 2ddT, a clo- 
sure condition is necessary, allowing one to construct the solution of 2ddT 
through an implicit system of algebraic equations, whose inversion is respon- 
sible for the wave breaking of an initial localized profile. In §3 we present 
the 1ST for the 2ddT equation (fT2l and use it to obtain the formal solution 
of the Cauchy for such equation. In §4 we obtain the longtime behaviour 
of the solutions of such a Cauchy problem, showing that the solutions break 
also in the longtime regime, and that such regime is essentially described 
by the longtime breaking formulae of the dKP solutions [26J , confirming the 
expected universal character of the dKP equation as prototype model in the 
description of the gradient catastrophe of two-dimensional waves. In §5 we 
characterize a class of RH spectral data allowing one to decouple and linearize 
the RH problem, generating a class of implicit solution of 2ddT parametrized 
by an arbitrary real function of one variable. 



2 Nonlinear RH dressing 

In this section we introduce the vector nonlinear RH problem enabling one 
to construct large classes of solutions of the Lax pair ([7]) and of the 2ddT 
equations ([2]) and (fT2|) . 

Proposition. Consider the following vector nonlinear RH problem 
^(X) = ^(X) + R j (^(X) + p 1 (X),^(X) + u 2 (X)), Aer, j = 1,2 (18) 

on an arbitrary closed contour T of the complex A plane, where R(s) = 
s 2 ), R 2 (si, s 2 )) t are given differentiable spectral data depending on 
the second argument s 2 through exp(is 2 ) and satisfying the constraint 

{^i(Si,S 2 ),^ 2 (si,S 2 )}( SliS2 ) = 1, , lg v 

TZj(si,s 2 ) := sj + R j (s 1 ,s 2 ), j = 1,2, 

with 

{f,g}(si,s 2 ) ■= f Sl g S2 - f S2 g Sl ; (20) 

where Vj, j — 1,2 are the explicit functions 



vi\ ( (CiA + C 2 A~> - t - Ci< 



v 2 j V i In A + if 



(21) 
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and £ + = (£i~,£2~) T and £ = (£ x , £ 2 ) T are the unknown vector solutions 
of the RH problem (TTSj) . analytic rispectively inside and outside the contour 
T and such that £ _ — > as A — > oo. Then, assuming that the above RH 
problem and its linearized form are uniquely solvable, we have the following 
results. 

1) If 

*£+(O)=0 4 , 

it follows that tF^ = ^ + z? are common eigenfunctions of Li >2 : Li7r ± = 
L2'JT ± = satisfying the relations 

tnf> }(A,t) = A(7if A ?r 2 t t - ?if f 7T 1 ± A ) = z (23) 

and the potentials reconstructed through (|22|) . solve the 2ddT equa- 

tion (ED. 

2) In addition, if the variables (1,(2 are specified as in (fTTT) . if the RH data 
satisfy the additional reality constraint 



K (K(s) )=s, VseC (24) 



and T is the unit circle, then the eigenfunctions satisfy the following symme- 
try relation: 

tT(A) =7f+(l/A) (25) 

and (f) G R. 

Remark 1 The RH problem (ITS]) can obviously be formulated directly in 
terms of the eigenfunctions ff as follows: 

7r+(A)=^(7rr(A),7r 2 -(A))=7r7(A) + J R,(7rr(A),vr 2 -(A)), A G T, j = 1,2, 

(26) 

with the normalization 

tt-(A) = z7(A) + 0(A" 1 ), |A| » 1. (27) 

— # 

Remark 2 The dependence of R on s 2 through exp(is 2 ) ensures that the 
solutions £ of the RH problem do not exhibit the In A singularity contained 
in the normalization (|21|) . It follows that the eigenfunctions contain the 
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In A singularity only as an additive singularity, while 7iy do not exhibit such 
singularity. 

Remark 3 Before adding the closure conditions (I2"2"j) . the solutions £ of the 
RH problem depend, via the normalization (I2TI) . on the undefined fields 0j 
and (f>£ 1} through the combination (t + Ci0Ci)> then the two closure condi- 
tions (|22|) must be viewed as a nonlinear system of two algebraic equations 
for (p t and 0^ defining implicitly the solution , (p t of the 2ddT equation. 
Therefore, as in the dKP case [26], we expect that this spectral features be 
responsible for the wave breaking of localized initial data evolving according 
to the nonlinear wave equation ([3]). Details on how two-dimensional waves 
evolving according to the 2ddT equation break will be presented elsewhere. 
An alternative closure, perhaps useful in the reality reduction case described 
in part 2) of the above Proposition, is given by the equations 



With this closure, indeed, we obtain a system of algebraic equations involving 
t + z<f) z , its imaginary part and 4> t . 

Remark 4 The symmetry relations ff25l) are a distinguished example of the 
following symmetry of the common eigenfunctions of the Lax pair (JTj), when 
Ci — z i C2 = z as in ( TTTj) and <fi G K: 

if ip(X) is a solution of |?|) ; then ipil/X) is a solution too. 

Remark 5 For part 2) of the above Proposition, when the contour V is 
the unit circle, the RH problem is characterized by the following system of 
nonlinear integral equations for £• (A), |A| = 1 (having parametrized A as in 



lm(t + z(f) z ), 



(28) 



2 



^2 + (0) = <Pt- 



my- 




(29) 



and the closure conditions (1221) read 




/ dde- iB R 2 ( (ze~ ie + ze i6 )v - t - z<f> z + £r(e" ie ), 



(30) 
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2vr / 

o V (31) 



Proof. For part 1), we first apply the operators Li 2 to the RH problem f l26lh 
obtaining the linearized RH problem Ljif + = ALjif" , where A is the Jacobian 
matrix of the transformation (|26|) : Aij = dlZi/dsj, i,j = 1,2. Since, due to 
the normalization (|27() . Lj-K~ — ► as A — ► oo, it follows that, by uniqueness, 
ff^ are common eigenfunctions of the vector fields L\^'- L\^^ = and, 
consequently, that <f> t , <p z are solutions of the 2ddT equation (J2]). Then the 
eigenfunctions exhibit the following asymptotics: 



l 2 

where 



tx- = \zv -t- z<j) z + X' 1 (zv + a"v~ x ) + 0(A~ 2 ), |A| >> 1, 

= X^zv -t -zfa + X(zv + a + v~ l ) + 0(A 2 ), |A| << 1, 
ti 2 = ilnX + - iX-^.v- 1 + 0(\- 2 ), \X\ » 1, 
tt+ = iln X - if + iA^ 1 + 0(A 2 ), |A| « 1, 



(32) 



<h = -4 = - (*e*)*- (33) 

implying the closure conditions fl22l . fl28l) . Since, from fl26l) . {n^ , vr^}( t)A ) = 
{^i,^2}7r-{7riT> 7r 2"}(t ) A)) A G T, equation ([191) implies that {^f ,^t}(t,\) = 
{fti, 7T2"}(t,A), A G T; i.e., the Poisson brackets of the ± eigenfunctions are 
analytic in the whole complex A plane. Since {^1,^2} -* i as A -> 00, it 
follows that {7r^,7r^} = {^1,^2} = For P ar ^ 2), applying 7£(-) to the 
complex conjugate of the RH problem (1261) and using the reality condition 



( 1241) it follows that vr + (A) = 7T~(A), |A| = 1. By the Schwartz reflection 
principle, it follows the symmetry (j25D and, using the Lax pair ([7]), the reality 
condition </> G R. □ 



3 The Cauchy problem for the 2ddT equation 

In this section we present the formal solution of the Cauchy problem for the 
wave form of the 2ddT equation: 

(e*) t = (ft xx + 4> yy , x, y G R, t>0, cf)(x, y, t) G R, 
(f>(x,y,0) = A(x,y), (f> t (x,y, 0) = B(x,y). 
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where the assigned initial conditions A(x, y), B(x, y) are localized in the (x, y) 
plane for x 2 + y 2 — > oo. To do it, we use the 1ST for vector fields developed 

in pa EQl Ell I2SJ. 

In this respect, we recall two basic facts: since the Lax pair of the 2ddT 
is made of vector fields, i) the space of eigenfunctions is a ring (if fi and f 2 
are eigenfunctions, any differentiable function F(fi, f 2 ) is an eigenfunction) ; 
ii) since the vector fields are also Hamiltonian, the space of eigenfunctions is 
also a Lie algebra, whose Lie bracket is the Poisson bracket (T^) (if f\ and 
fi are eigenfunctions, also /^(a,*) is an eigenfunction). 

Multiplying the first and second equations of (J7|) (with (j = z, (2 = z as 
in (ITT]) ) by A -1 and A respectively, then adding and subtracting the resulting 
equations, one obtains the equivalent and more convenient Lax pair: 

CnJ> := A</> 2 - - A"% - ( - 2v t + A^f + A-^A^a = 0, (35) 

:= Vt + ^ (A^- + A"V*) - t 1 (A0,-t - A- 1 ^*) A^a = 0, (36) 

where the first equation must be viewed as the spectral problem (in which v t 
shall be replaced, in the direct problem, by (4> zz -/2)exp(—(fi t /2), due to (TT2T) ) 
and the second equation as t-evolution of the eigenfunction. 

Eigenfunctions and spectral data. Now we introduce the Jost and an- 
alytic eigenfunctions for the spectral problem (|35|) . Since the associated 
undressed operator: Xd z — \~ 1 d z coincides with the undressed operator of 
the spectral problem for the (2 + l)-dimensional self-dual Yang-Mills equa- 
tion [44.J, the construction of the Jost and analytic Green's functions is taken 
from there. 

We define Jost eigenfunctions of the spectral problem fl35l) on the unit 
circle of the complex A plane, using the parametrization (|17[) . Introducing 
the convenient real variables £, 77, 9' as follows: 

£ = cos 9 x + sin 9 y, 

r] = — sin 9 x + cos 9 y, (37) 
9' = 9, 

the Lax pair f[3"5l) . f[3"6"l) becomes 

Clip := Tp v - \ [- (0 CC + <p m ) v- 1 + <p & ] (77^ - £tp v + Vfl') = °> ( 38 ) 

t 2 ^ : = ^ + v- 1 ^ - (v- l ) ri ( V ^ - + fa) = 0. (39) 



9 



A convenient basis of Jost eigenfunctions are the solutions f\ and /2 of equa- 
tion (1381) satisfying the boundary conditions 

f(£,r),e') := ( f ^ V J 9 2 ) - ( 1 ), as 77 - -oo; (40) 



they are characterized by the linear integral equation 

/=(!') + 5 / <V + ^vK 1 + M {ih - Urf + fo')- (4i) 

It follows that rj, 9) and /2(C) Vi ®) ~ ® are 27r-periodic in 9. 

The 77 — > 00 limit of / defines the scattering vector <?(£, 6>) = #), cr 2 (£, #)) T 

as follows 

f(£,ri,0)^£(Z,0)= ( l)+*(£,0), as 77 - 00; (42) 



namely: 



a(£, 9) = \\di ] + ^JiT 1 + <^J (77/^ - + £). (43) 



Also the scattering vector is 27T-periodic in 9: cr(£, 6 + 2n) = <?(£, 9); i.e., its 
dependence on the second argument 9 is through exp(i9). 

The analytic eigenfunctions of ( |35l) are defined instead via the integral 
equations: 

I j c dz' A dz'G^z -z'^-z', A) + A%i + A" 1 ^] At/^', A), 

(44) 

where G are the analytic Green's functions 

G ± (z,zA) = T- n sgn(l - |A|) = ±1 (45) 

such that AGf — X^Gf = 5(z), reducing, on the unit circle |A| = 1, to 

G ± (z,z,X) = T -—t— , 0<e«l. (46) 

7T ^ =F 7677 
10 



Since G + and G~ are analytic respectively inside and outside the unit 
circle of the complex A plane, then (ipf, ip^) and (V'iTj ^2") are a ^ so analytic, 
respectively, inside and outside the unit circle of the complex A plane, af- 
ter subtracting their singular parts, given respectively by (A _1 £i>, z In A) and 
(Azt>,ilnA), as it can be seen by solving the integral equations (J4"4"j) by iter- 
ation or from the following A - asymptotics: 



V>r = Xzv - z(f) z + A" 1 (zv + a~v~ l ) + 0(A~ 2 ), |A| >> 1, 
ipf = X~ x zv -z(f) z + X (zv + a + v- 1 ) + 0(A 2 ), |A| << 1, 
$2 = iln A + if - iA~ + 0(A- 2 ), |A| » 1, 
= ilnX - + iX^zV- 1 + 0{X 2 ), \X\ « 1, 

where a are defined in fl33l) . 

In addition, equations (1460 imply the limits 



(47) 



G+(z- z',z- z', X) -»• -J|=^±fc, as 77^^00, 
- z', z - z', A) -» ; f _ f L , as 77 ^^00. 



(48) 



Therefore, on the unit circle |A| = 1, the 77-^—00 limit of (ipt^t) and 
(ipijipz) are analytic respectively in the upper and lower parts of the com- 
plex £ plane, while the 77 — > 00 limit of {ipt-i^t) and (ipi,^) are analytic 
respectively in the lower and upper parts of the complex £ plane. This mech- 
anism, first observed in [H], plays an important role in the 1ST for vector 
fields (see [IS], [2D], [21], [25]). 

Since the Jost eigenfunctions / = (/1, /2) 71 are a good basis in the space 
of eigenfunctions of the spectral problem fl38|) for |A| = 1, one can express 
the analytic eigenfunctions in terms of them through the following formulae, 
valid for |A| = 1: 

k = ^ k (f)=f+X k (fi,h), (49) 

defining the spectral data as different iable functions of two arguments. 
In the 77 — > — 00 limit, equations (|49l reduce to 



ii.m _'0 ± - ( I ) = x k (i.»). (W) 



17- 



implying that i) x + (£,#) and 0) are analytic in the first variable £ 

respectively in the upper and lower half parts of the complex £ plane, and ii) 
X^iCi &) are 27r-periodic in 6: X^i^i 9 + 2vr) = ^(^ 6) (their dependence on 
the second argument 9 is through exp(i6)). 
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At rj — > oo, equations reduce to 

lim^- ( I) =c?+x ± (e + a 1 ,^ + a 2 ). (51) 



2t 



Applying the operator Ld£ J |i e - J ( w €+ n6 '). ) n G Z to equations flBTl) and 

o 

using the above established analiticity properties in £ and the 2-7r-periodicity 
in we obtain the following linear integral equations connecting the (Fourier 
transforms of the) scattering data a to the (Fourier transforms of the) spectral 
data x ± - 



X±(u;, n) + H(±u) I a(u, n) + / eL/ £ n')Q(u', n', u, n) ) = 0, 

V K n'=-oo / 

(52) 

where if is the Heaviside step function and 

Q(u/,n',u,n) = / § Jg e ^'-")+K-«W ( e ^(€.«)+nW€,e)) _ i) ? 

R 

^(w.n) =/de/ £ e-^+^x + (e^), (53) 

R 



— » 

At last, eliminating, from equations ( H9l) . the Jost eigenf unctions /, one 
obtains, through algebraic manipulation, the following vector nonlinear RH 
problem on the unit circle of the complex A plane: 

$x =Ki(fo,fc) =ipi + Ri(ipi,ip2), |A| = 1, ( v 

^+ = ^ 2 (^-,^-) = ^-+i? 2 (^-,^-). 10 ' 

We remark that the 27r-periodicity properties of the scattering data x^(£, 0) 
in the variable imply that the dependence of R on the second argument 
S2 is also through exp(is2), to guaranty that the In A singularity is just an 
additive one for ip 2 > an d is absent for ^ . 

Recapitulating, in the direct problem, at t — 0, we go from the initial 
conditions 4>, 4>t of the 2ddT equation to the initial scattering vector cr(£, 0); 
from it we construct, through the linear integral equations ( l52l) . the scat- 
tering data %(£,9) and, through algebraic manipulation, the RH spectral 
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data R(s) = (R%(si, s 2 ), ^2(^1, s 2 )). In the inverse problem, one gives the 
RH spectral data R(s) and reconstructs the vector solutions ip^ of the RH 
problem f|54|) . defined by the normalization: 

t- = ( {zx Y^r^t" z4 " ) + °" (a " ) ' iai >y 1 <55) 

At last, the closure conditions 



lim \(iip 2 +hiA) = 4> z e 2 . 

A^oo 

lim (iibt + In A) = — , 



(56) 



consequences of the asymptotics ()47j) . allow one to reconstruct the solution of 
the 2ddT equation through the solution of a system of two algebraic equations 
for (f) t and <j) z . 

Time evolution of the spectral data. To construct the t-evolution of 
the spectral data we observe that / and ip^, eigenfunctions of the spectral 
problem fl35l) : L\j = = 0, are solutions of the following equations 

involving the second Lax operator £ 2 / — — (1, 0) T , implying the 

following elementary time evolutions of the data: 

*)=<?(£-*, 0,0), x ± (^9,t) = x ± (^-t,e,o), 

In addition, it follows that the common Jost eigenfunctions J and the com- 
mon analytic eigenfunctions tF of the Lax pair (I35l) . (l36"l) are obtained from 
/ and simply as follows: 

J:=/-t(l,0f 
tt ± := ^ - t(l,0) 



± — ,T± _ in ml ( 58 ) 



— * 

It is easy to verify that the analytic eigenfunctions 7? , the RH data R(s) 
and the associated RH problem of this section coincide with those appearing 
in the dressing construction of §2. 

Hamiltonian constraints on the data. The Hamiltonian character of the 
2ddT dynamics implies the following formulae for the Poisson brackets of the 
relevant eigenfunctions: 

Ul, ^2>(A,t) = {Tlf, 7r 2 = }(A,t) = i (59) 
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which, in turn, imply that the transformations s — > ^C ± (s) and s — > lZ(s) are 
canonical: 

{/Cf , Kf} (suS2) = {K u K 2 } (suS2) = 1. (60) 

To prove fl59|) . one first shows that J 3 := {Ji, J2}{x,t) — ► i as 77 — > — oo, 
^3 := {^r. ^2 }(A,t) — > * as A — > oo, 7r^ := {vr^, tt^ -> i as A -> 0. Since 
the vector fields are Hamiltonian, J 3 , 7T 3 are also common eigenfunctions, 
and equations (1591) hold, by uniqueness. Equations (JB"D]) are consequences of 
and of the relations 



7T 



± 



(J) , 7T+ = 7? (tT) . (61) 



Reality constraints. The definition (158b ) and the condition e K imply 
the symmetry relations 

7T(A) = 7f+(l/A); (62) 

consequently, from ( fT8l) . the reality constraint ( 1241) on the RH data holds 
true. 

Small field limit and Radon Transform. As for the 1ST of the heavenly 
[2"U] and dKP [23] equations, in the small field limit |0|, \4> t \ « 1, the direct 
and inverse spectral transforms presented in this section reduce to the direct 
and inverse Radon transform [35]. Indeed, the mapping from the initial data 
{A(x, y), B(x, y)} to the scattering vector a reduces to the direct Radon 
transform: 

0) ~ I / ( i ) [ - ( 9 I + 9 v) A «C, V, o)MZ, v, o))+ 
dsB{x((,7i,0),y((,7i,0)j\dri, 

x (£> Vi @) = £ cos — r/sm 9, y(£, rj,9) = £ sin 9 + r\ cos 9, 



(63) 



while the spectral data and R are constructed from a as follows: 

^(e, ^) ~ -pf e), ~ -in^ic, 9), (64) 

where and 7l!g are rispectively the (±) analyticity projectors and the 
Hilbert transform in the variable £: 

= ^ /„ e^5) 9K ''' H{9(?) := ^ p / i^ 9ia (65) 
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At last, the first of the closure conditions (155]) of the inverse problem reduces 
to the inverse Radon transform 



2?r 2tt 

m P r «MF^)" 2 ^ (66) 



o 

£ = x cos + ysin^, 
that can be shown to be equivalent to the well-known Poisson formula 

<j>(x, y, t) = d t j R2 d 4^L[x -x',y- y', t)A(x>, y>) + / Ra d 4£fh{x - x', 
y-y',t)B(x',y'), 

(67) 

where 

L(x,y,t):=^=f£=z£± (68) 
y/t z — x z — y z 

and H(-) is the Heaviside step function, describing the solution of the Cauchy 
problem 

4>tt = 4>xx + <P yy , x,y eR, t>0, (f>(x,y,t) el, , gg v 
(f>(x,y,0) = A(x,y), (f) t {x,y, 0) = B(x,y). 

for the linear wave equation in 2+1 dimensions. 



4 The longtime behaviour of the solutions 

In this section we show, as it was done in the dKP case [26], that the spectral 
mechanism causing the breaking of a localized initial condition evolving ac- 
cording to the 2ddT equation is present also in the longtime regime. We will 
actually show that the longtime breaking of the 2ddT solutions is essentially 
described by the longtime breaking formulae of the dKP solutions found in 
|26j ; this is an important confirmation of the expected universal character 
of the dKP equation as prototype model in the description of the gradient 
catastrophe of two-dimensional waves. 

We remark that it is clearly meaningful to study the longtime behaviour 
of the solutions of the 2ddT equation only if no breaking takes place before, 
at finite time. In this section we assume that the initial condition be small, 
then the nonlinearity becomes important only in the longtime regime and no 
breaking takes place before. 



15 



Motivated by the longtime behaviour of the solutions of the linear wave 
equation u tt = u xx + u yy , localized, with amplitude 0(t~~2), in the region 
a/x 2 + y 2 — t = 0(1), we study the longtime behaviour of the solutions of 
the 2ddT equation in the space-time region 

z =*±i e ia , a,reR, a = 0(1), t » 1, (70) 

implying that 

r = yjx 2 + y 2 — t, a = arctan — . (71) 

x 

Substituting (17D1 into the integral equations (129|) and keeping in mind that, 
in the longtime regime, <p t is small, so that, f. L, v ~ 1 + 0t/2 + 2 /8, we 
obtain 



(a) - £ / ^4%'-*) ^ ( - 2t sin2 + r cos ^' - «)+ 

^ cos(fl' - a)&(l + f ) - # 2 + £r(e _ *'), 0' + frCe"*')) ~ 0, j = 1, 2. 

(72) 

Since the main contribution to these integrals occurs when sm((6'—a)/2) ~ 0, 
we make the change of variable 6' = a — /// ^/i, obtaining 



R l-(l=Fe)e V V / (73) 



a + $\e< a - J *>\ )~0, J = 1,2, 



where 



X := r + + |# - # 2 . (74) 



If \9 - a| » r 1 / 2 , equations (O imply that £f{\) = O^ 1 / 2 ): 



(A) r? 7^-r-7 m\ f <^'#i f-^+^ + ^r(e" 



« + e 2 -(e- 4(Q -^l ), j = 1,2. 



(75) 



If, instead, - a = -/i^ 1 / 2 , M = 0(1), then ff (A) = 0(1): 



IR \ 



« + e 2 -(e- l(Q -^l ), j = 1,2. 



(76) 
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Therefore it is not possible to neglect, in the above integral equations, £~, j = 
1, 2 in the arguments of Rj, j = 1, 2; it follows that these integral equations 
remain nonlinear even in the longtime regime. 

At last, using equations (1751) . the asymptotic form of the closure condi- 
tions read, for t » 1: 



« + £ 2 e ^ 



(77) 



2 = _ e — ^(1 + f ) (1 + ©(r 1 )) . (78) 

Comparing (jTTjl and (1751 . and using flTOj) . we infer that z0 2 ~ — ^^ t (l + ^). 
Using this asymptotic relation in (ITT)) , we finally obtain the following result. 

In the space-time region 

2 =*±V a , a,r6l, t»l, 

X : =r+ (t + r)<p t + f = (79) 

V^+F - * + v^+rt + f 4>i = o(i), 

the longtime t » 1 behaviour of the solutions of the 2ddT equation {exp4>t)t = 
(fixx + (fiyy is described by the following implicit (scalar) equation: 

<f>t = ^ t F ( v^+F - t + y/rf + lfh + f 0?, arctan f ) + o (^) , (80) 
where F is given by 



F(X,a) = ~^-J dfi'R 2 (- £ + X + a 1 (//;X,a),a + a 2 (/i , ;X, 

R V 

and aj(/x; X, a), j = 1,2 are the solutions of the integral equations 

aj {n; X,a) = £rf j^z^Rj + a^'; X, a), a+ 

a 2 (//;X,a<)Y j = 1,2. 



a. 



(82) 



Outside the asymptotic region (179")) the solution decays faster. 

We first remark that, since <p t = 0(t~ 5 ), the condition X = 0(1) implies 
that r = \Jx 2 + y 2 — t = 0(y/i); it follows that, in the longtime regime t » 
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1, the solution of the Cauchy problem for the 2ddT equation is concentrated, 
with amplitude 0{t~^), in the asymptotic region y/x 2 + y 2 -t = 0(y/i). We 
also remark that the asymptotic solution (180]) - (182]) is connected to the initial 
conditions of the Chauchy problem through the direct problem presented in 
the previous section. 



5 A distinguished class of implicit solutions 

In this section, in analogy with the results of [26], [28], we construct a class of 
explicit solutions of the vector nonlinear RH problem ffl8l) and, correspond- 
ingly, a class of implicit solutions of the 2ddT equation parametrized by an 
arbitrary real spectral function of one variable. 

Suppose that the two components of the RH spectral data R in (fl8|) are 
given by: 

R J (s 1 ,s 2 ) = (-iy + Hf(e s ^), j = 1,2, (83) 

in terms of the single real spectral function / of a single argument, depending 
on si and s 2 only through their sum. 
Then the RH problem (126]) becomes 

and the following properties hold. 

i) The reality and Hamiltonian constraints (I24"j) and (fl9l) are satisfied. 

ii) irf + = + n 2 . Consequently, using the analyticity properties of the 
eigenfunctions, it follows that the functions A + and A~, defined by 

A ± :=7if + vif -(zA + zA-^u-ilnA, (85) 

are analytic respectively inside and outside the unit circle of the A-plane and 
satisfy the equation A + = A~; therefore they are equal to a constant in A. 
Evaluating such a constant at A = and at A — > oo, we obtain the following 
equalities 

A+ = A" = -t-zfe-i— = -t-z<p z + i— . (86) 

2 2 
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This implies that 

i) the solutions of the 2ddT equation generated by the above RH problem 
satisfy the linear (2 + l)-dimensional PDE 

(fit = i(z<p s - zcfi z ) (87) 

and, substituting in (fl2"|) the expression of (fit in terms of <f> z , (fig given in (!H7|) . 
one obtains the following nonlinear two dimensional constraint: 

i{zd- z - zd z ) ( e <(^«-^»)) = <j> xg (88) 

on the solutions of 2ddT constructed by the above RH problem. 

ii) rti + = Tti + is the following explicit and elementary function of A: 

w(X) := + 7^2 = + n 2 = (zX + zX~ 1 )e~^ + i In A — t — zcfi z + i^. 

(89) 

hi) Since, from (1891) . tt 1 + 7r 2 = w(X) is an explicit function of A, the vector 
nonlinear RH problem (1841) decouples into two scalar, linear RH problems: 

7T+ = 7T7 + if (e«W) , 

- 2 + = vr 2 - - i/ ( e «W) , (90) 
whose explicit solutions are given by 

ef(A) = v- (1 Se^ / (^ (V) ) . J = L2, (91) 

where = — Uj, and the closure conditions (122]) read 

^ = ~<f dXf(e^), cfi t = ±-l ^/(e«W). (92) 
2m J| A | =1 2m yj A | =1 A 

Although the RH problem ( 1901 is linear, since w(X) in (l89l) depends on the 
unknowns <fit,<fi z , the closure conditions (1921) are a nonlinear algebraic system 
of two equations for the two unknowns <fit,<fi z , defining implicitly a class of 
solutions of the 2ddT equation parametrized by the arbitrary real spectral 
function /(•) of a single variable. 
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